Inseparability criteria for bipartite quantum states 
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We provide necessary and sufficient conditions for the partial transposition of bipartite harmonic 
quantum states to be nonnegative. The conditions are formulated as an infinite series of inequalities 
for the moments of the state under study. The violation of any inequality of this series is a sufficient 
condition for entanglement. Previously known entanglement conditions are shown to be special 
cases of our approach. 
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Entanglement plays a key role in the rapidly devel- 
oping field of quantum information processing. In this 
context it is important to provide methods for character- 
izing entangled quantum states on the basis of observable 
quantities. However, already in seemingly simple cases 
this problem turns out to be rather complex. Even for a 
two-party harmonic-oscillator system so far there exists 
no complete characterization of entanglement to be used 
in experiments. 

For characterizing entanglement, that is inseparability, 
of the density operator of a bipartite continuous variable 
system, one may use the Peres-Horodecki condition 0, 
I2, y ■ A sufficient condition for entanglement consists in 
the negativity of the partial transposition (NPT) of the 
quantum state of the two-party system. To characterize 
NPT for such a system completely, however, to our best 
knowledge is still an unsolved problem. 

A sufficient condition for the NPT has been proposed 
by Simon |3[ . It is based on second-order moments of po- 
sition and momentum operators. For the special case of 
Gaussian states the resulting entanglement criterion has 
been shown to be necessary and sufficient. Another in- 
separability condition based on second moments has been 
derived without explicitely using the NPT condition [^. 
This condition is also complete for the characterization 
of entanglement of two-mode Gaussian states. The lat- 
ter approach has been extended to special higher-order 
moments ill] and even to more general operator func- 
tions 

The complexity of the problem under study may be- 
come clear when we go back to a related but simpler prob- 
lem. It consists in the characterization of nonclassical 
effects based on the negativity of the Glauber-Sudarshan 
P-function. Only recently the problem was solved of how 
to characterize the nonclassicality in terms of observable 
quantities. This requires an infinite hierarchy of condi- 
tions formulated either in terms of characteristic func- 
tions 0] or in terms of observable moments 0, Q| . 

In the present contribution we will further develop the 
concept of the complete characterization of single-mode 
nonclassicality with the aim to characterize the entan- 
glement of bipartite continuous variable quantum states. 
Based on the NPT condition we derive a hierarchy of 



necessary and sufficient conditions for the NPT in terms 
of observable moments. Even though this only leads to 
sufficient conditions for entanglement, it can be applied 
to a variety of quantum states. It does not only con- 
tain as a special case the condition of Simon 0| , but also 
other types of inequalities ^ tOh ilt the sake of 

strictness of the conclusions we restrict our attention to 
a two- mode harmonic oscillator. Hence the derived con- 
ditions are useful for radiation modes, motional states 
of trapped ions and related systems. The further exten- 
sion of our approach to more general continuous variable 
systems seems feasible, but this requires some further in- 
vestigations. 

A bipartite quantum state g is called separable if it is 
a convex combination of factorizable states. 
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where and pj"'' ^-re states of the first and the second 
part respectively. The p„ are nonnegative numbers sat- 
isfying J2n=aPn = 1 ^^"^ ^^'^ series in is convergent 
in the trace norm. A state which cannot be represented 
in the form is called entangled. 

It is well known that the full transposition of any (one- 
partite or multi-partite) state is again a quantum state. 
But in the multipartite case, along with the notion of 
the full transposition, there is a notion of partial trans- 
position, which transposes only some parts of the state 
and leaves the others unchanged. In contrast to the full 
transposition, partial transposition does not necessarily 
map a quantum state onto a quantum state. Clearly, 
the partially transposed density operator g^"^ of a sep- 
arable state g, Eq. is again a quantum state (we 
consider transposition of the second part; our results are 
unchanged if we transpose the first part): 
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smce g2 is a quantum state for all n and so is g^"^ ac- 
cording to the expansion (j^J . This property of separable 
states allows us to use the NPT as a sufiicient condi- 
tion for inseparability. This condition is referred to as 
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Peres-Horodecki condition 0, 13, • In general the par- 
tial transposition of a density operator may no longer be 
nonnegative defined. It is this property that we will use 
for the derivation of our hierarchy of conditions. 

Remember that a Hermitian operator A is called non- 
negative if 



(3) 



holds for all states The left hand side of this inequal- 
ity can be written as 



= Tr(i|^}(V'| 



(4) 



The nonnegative operator \ip){ip\ can be represented in 
the form / where / may be written as 



/ = |vac)(V'|. 



(5) 



The state |vac) is the vacuum state, in our bipartite case 
it reads, in the Fock basis, as |vac) — |00). Any pure 
bipartite state \ip) can be expressed as 



l^> =5^100), 



(6) 



for an appropriate operator function g = ^(d, 6), where d 
and b are the annihilation operators of the first and the 
second mode respectively. Hence the operator (O is of 
the form 



/= |00)(00|.g. 



(7) 



The two-mode vacuum density operator may be ex- 
pressed in the normally-ordered form. 



|00)(00| = : exp( 



-d^d 



StS) 



(8) 



with : • • • : denoting the normally-ordering prescription. 
By combining Eqs (O and ^ it immediately follows 
that the normally-ordered form of the resulting opera- 
tor / exists. Hence we conclude: a Hermitian operator 
is nonnegative if and only if for any operator f whose 
normally- ordered form exists the inequality 



{Pf)A = Tr(i/t/) ^ 



(9) 



is satisfied. 

Now we can use the inequality © to apply the Peres- 
Horodecki condition, which gives the following result. 
For any separable state g the inequality 

(/t/)PT ^ (/t/).,, = Tr(^PT/t/) ^ (10) 

holds true for any operator f whose normally- ordered 
form exists. According to the latter assumption the op- 
erator / can be expanded as 
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Upon substituting this expansion into the inequality 1101) 
we get the condition 
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^pqrs^i^^'^kl-^pqrs.nrakl 0, (1^) 



for the moments of the partial transposition 

M, 



(13) 



pqrs^n^kl = (dt?dPdt"d"6t«fo'-fotfe5/)PT, 

The left hand side of the inequality H12|l is a quadratic 
form with respect to the coefficients Cnmki of the expan- 
sion (|ll(l . The Silvester criterion states that the inequal- 
ity (|12|l holds for all Cnmki if and only if all the main 
minors of the from (|12|l are nonnegative. To write these 
minors explicitly we need to find the relation between the 
moments H13(l of the partially transposed state and those 
of the original one. 

It is convenient to number multi-indices {n,m, k,l) by 
a single number so that the moments Mpqrs,nmki be num- 
bered by two indices like matrix elements are. The prob- 
lem of numbering the multi-indices (n, m, fc, I) is equiv- 
alent to the ordering of the set {(n, m, k, I)} of all such 
multi- indices. Then we can use the ordinal number of a 
multi-index (n, m, k, I) in the ordered sequence of multi- 
indices. We use the following order: for any two multi- 
indices u = (n, m, fc, I) and v = (p, q, r, s) 



u <v 4^ 



|u| < |v|, or 

lul = Ivl and u <' v. 



(14) 



where |u| — n -\- m + k + I and u <' v means that the 
first nonzero difference i — k,s — l,p^n,q — m is positive. 
The resulting ordered sequence of the moments starts as 
follows: 



1, (d), (dt), (S), (6t), (d2>, (dtd), (dt2), (d6), (dtfo), 
{P),{ab^),{a^b^)Ab^b),{b^^),.... 



(15) 



One can easily check that the moments Mpqrs,nmkh 
Eq. (|13|l . can be expressed in terms of the moments of 
the original state as follows 

{a^iaPa^''aH^'b'bHY^ = {a^-JaPa^''dJ^bHHH'). 

Note that these moments are linear combinations of 
normally-ordered ones, {a^'^a^b^^U) , according to the 
following formula for the antinormally-ordered product 



n\m\ 



inin(n,m) 

d"d^™ = 

^ k\{n - k)\{m - k)\ 



dt™"'^d"^^ (16) 



The nonnegativity of all the minors of the quadratic 
form (|12|l can now be written as 
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where the moments Mij defined in Eq. (|13|l are now given 
as 



(18) 



with (n, m, fc, I) and (p, g, r, s) being the ith and the jth 
indices in the ordered sequence \ib\ . Now we can exphc- 
itly formulate a necessary and sufficient condition for the 
positivity of partial transposition. The partial transposi- 
tion of a bipartite quantum state is nonnegative, if and 
only if all the determinants 



Dn = 



1 
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(atfl) 
(a2) 

(ab) 
(aSt) 



(at) 
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(19) 



are nonnegative, that is 

\fN : Dn> 0. 



(20) 



The other way around, if there exists a negative deter- 
minant, 



3N : Dn <0 



(21) 



the NPT has been demonstrated. Hence the condi- 
tion \21]) is a necessary and sufficient condition for the 
NPT. That is, it provides a sufficient condition for the 
state under consideration to be entangled. This result 
is the central finding of the present contribution. It is 
important to stress that all the moments occurring in 
our conditions can be determined experimentally. This 
requires a straightforward extension of the method, as 
proposed for measuring the moments of a single-mode 
radiation field 0|, to the case of two modes. 

Even though the above conditions filfl are necessary 
and sufficient for NPT, for practical applications it may 
be useful to derive other types of conditions. They may 
be more efficient for verifying entanglement of particu- 
lar quantum states in practice, violations of these new 
conditions may be found easier than the negativities of 
the determinants H19|l. For this purpose we consider the 
following two possibilities: 

(i) Restricting the generality of the operator l|ll|l by 
restricting the generality of the coefficients Cnmki 
leads to other conditions for inseparability. For ex- 
ample, setting sonic coefficients Cnmki to zero leads 
to conditions for the nonnegativity of subdetermi- 
nants of _Djv, which are obtained by canceling rows 
and columns with the same indices. 

(ii) Identifying some coefficients Cnmki leads to another 
way to formulate new conditions for NPT by re- 
ducing the number of independent variables of the 
form H12|l . For example, the quadratic form derived 



for the operator / = ci -1-020-1-036 depends on three 
independent (complex) variables ci, C2 and C3. If 
we identify the last two coefficients C2 = C3, i.e. 
f — ci -\- C2{a -\- b) , we obtain a quadratic form of 
two variables ci and C2. The resulting sufficient 
conditions for NPT can be formulated in terms of 
the nonnegativity of determinants whose entries are 
combinations of the moments (|18|l . 

Examples for the usefulness of both types of new condi- 
tions will be considered in the following. 

Let us compare our approach with some previously 
known results. We start with the separability condition 
found by Simon ^J- Define matrices Ai, A2, C and J by 
the expressions 



A, = 



((Ax,)2) ({Ai-„AK}) 
({Ax.,ApJ) ((Ap,)') 

/(AX1AX2) (AxiApa) 
(^(ApiA^a) (ApiAp2) 



(22) 



i*]^ Q ) . The Simon condition can be formu- 



and J 

lated as follows: for any separable quantum state the 
quantity S, 



(23) 



S = det Ai det A2 + Q + dot C 
Tr(yli JCJAaJC^j) - i (dot Ai det A2) ^ 0, 



is nonnegative. Since the position and momentum oper- 
ators Xi and Pi, i = 1,2, are linear combinations of the 
creation and annihilation operators, the condition 1231) 
can also be formulated in terms of the moments M,j , 
Eq. p8|l . By using our concept the left hand side of the 
inequality H23|l reads much simpler in terms of the mo- 
ments Mij. Simon's quantity S is nothing but our fifth 
order determinant D5: 



S = D. 



(24) 



Therefore the Simon condition (|23|) is exactly one condi- 
tion of our hierarchy (|2(J|I . 

Another separability condition has been formulated 
in 5j. For any non-zero real parameter r one may de- 
fine operators u and i) via 

u = |r|i;i -f r~-^i;2, v — \r\pi — r~'^p2. (25) 

Then, for any separable state the inequality 

((Au)2) + ((At))2) - (r2 + r-2) ^ (26) 

holds true for all r ^ 0. The condition H26|l can be simpli- 
fied as follows. We express the operators u and v as linear 
combinations of the creation and annihilation operators 
and substitute them into the inequality H2t)|l . Now its left 
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hand side can be minimized with respect to r. Then the 
inequality (|26() is found to be equivalent to 



Explicitly it reads as 



(AatAa)(A6^A6) ^ Re^(AaAS). 



(27) 



The inequality (|27() is related to a subdeterminant of ZJjv, 
to the one resulting from the operator / — ci + C2a + C36: 



1 (a) (St) 
(at) (flta) (atSt) 

(6) (aS) (St6) 



(28) 



The condition d ^ explicitly reads as 

(AatAa)(A6tA6) ^ |(AaAS)|^ (29) 

from which it becomes clear that the condition (|26|) is a 
weaker form of the nonnegativity of our determinant (|28|l . 

There have been various other results in the litera- 
ture. We can only briefly comment on some of them. A 
generalization of the inequality H26|) was obtained in 
The main result is not explicitly based on the nonneg- 
ativity of partial transposition, but in fact it is a con- 
sequence of it. Using our theory we can show that the 
main inequality follows from Eq. I|10(l for the operator 
/ = ci_-|- 02^1 -l- C3B1 + C4A2 + C5-B2 (in the notation of 
work which is readily obtained by our approach (ii). 
Another inequality has been obtained in It differs 
from the condition H26(l and its generalization ^ in that 
it gives a lower bound for the product of the variances 
of the operators (|25|l rather than for their sum. Though 
this inequality does not explicitly use NPT, in fact it is 
also a consequence of it. The condition provided in 
easily follows from the nonnegativity of d, Eq. H28|l. to- 
gether with the determinant corresponding to the oper- 
ator / = ci + C2{d + b) + C3(at -I- 6t). Yet another con- 
dition has been recently obtained in It is easy to 
verify that this condition is a weaker form of the nonneg- 
ativity of the determinant corresponding to the operator 
/ = ci + C2ab + cafltfot. 

Let us illustrate our theory for the example of an en- 
tangled quantum composed of two coherent states, 

\^-) = Af-{a, /3)(|«, /?)-!-«, -/?)) , (30) 

where the normalization is 

-1/2 



AA_(a,/3) = (2(l-e-2(l"l'+l/3|'))y 



(31) 



The violations of the conditions H23(l and l|26|l fail to 
demonstrate the entanglement of this state. However, 
we may apply the following simple subdeterminant: 



1 



(St) 
(StS) 



(a6t) 
(ab^b) 



(atfo) (fltfotfe) (fltafetfo) 



(32) 



3th(|c 



sinh (lap 



(33) 



which is clearly negative for all non-zero parameters and 
hence the entanglement is verified. 

Finally, we will consider another special kind of sepa- 
rability condition obtained by our approach (i). Let us 
consider the simplest subdeterminants of the Dn, namely 
those of second order. They correspond to the condi- 
tion (|10|l with / being a linear combination of two terms. 



/ = ciat"a"6tfe5i + c2d'iPa%H''. (34) 
Explicitely the resulting separability condition reads as 

(at^a^atn^m^ti^fe^tfc^/) {a^idPa^Pd%^'b'^b^''b') ^ 

. . . . 2 (35) 
(at"o"atPa95t«^,r^tfe^,i) . 

This leads to rather simple conditions including moments 
of higher orders. Note that the conditions H35|l generalize 
those derived recently in • 

In conclusion, we have derived a hierarchy of condi- 
tions that is necessary and sufficient for the negativity 
of the partial transposition of the density operator of a 
bipartite harmonic oscillator. In fact this provides a suf- 
ficient condition for entanglement in terms of observable 
quantities, that can be used for a variety of practical ap- 
plications. We have also demonstrated that several pre- 
viously known conditions appear to be special cases of 
our approach. 
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